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Abstract 

An equilibrium of a planar, piecewise-C 1 , continuous system of differ- 
ential equations that crosses a curve of discontinuity of the Jacobian 
of its vector field can undergo a number of discontinuous or border- 
crossing bifurcations. Here we prove that if the eigenvalues of the Ja- 
cobian limit to A^iiu^ on one side of the discontinuity and — A^iio^ 
on the other, with \l, \r > 0, and the quantity A = Xl/^l — Xr/u>r 
is nonzero, then a periodic orbit is created or destroyed as the equi- 
librium crosses the discontinuity. This bifurcation is analogous to the 
classical Andronov-Hopf bifurcation, and is supercritical if A < and 
subcritical if A > 0. 

PACS: 02.30.Oz; 05.45.-a 



1 Introduction 

Many dynamical systems can be modeled by a system of differential equations 
z = F(z) with vector field F : M — ► MJ 1 on an n-dimensional manifold M. 
Such a system is piecewise- smooth continuous (PWSC) if F is everywhere 
continuous and is smooth except on the boundaries of countably many regions 
where it has a discontinuous Jacobian, A = DF. These boundaries are called 
the switching manifolds. 



* We acknowledge support from the National Science Foundation through grant DMS- 
0202032. 
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Piecewise-smooth dynamical systems are encountered in a wide variety of 
fields. Examples include vibro-impacting systems and systems with friction 
[Bro99l IWDKOOj . switching circuits in power electronics [BVOlj . relay control 
systems [ZM03] and physiological models [KS01] . 

Many bifurcations caused by discontinuities have been previously studied. 
Feigin first studied period-doubling bifurcations in piecewise continuous sys- 
tems and gave them the name "C-bifurcations" for the Russian word svejnye 
for "sewing," and this term is often used more generally for any bifurcations 
that result from a discontinuity |DBFHH99l IDBGIV021 ILN041 ILei06j . The 
bifurcations caused by the collision of a fixed point of a mapping with a 
discontinuity were studied by |NY92| and given the name "border-collision" 
bifurcations. Other examples of C-bifurcations include discontinuous saddle- 
node bifurcations (see $2]), grazing bifurcations, where a periodic orbit touches 
a switching manifold, and sliding bifurcations, where a trajectory moves for 
some time along the switching manifold. 

In this paper we study the bifurcation analogous to the Andronov-Hopf 
bifurcation that arises when an equilibrium encounters a switching manifold. 
Recall that in the classical Hopf bifurcation, a periodic orbit is generically 
created when an equilibrium has a pair of complex eigenvalues cross the imag- 
inary axis |Kuz04t IMM76] . We will obtain a similar result for the case that 
an equilibrium of a planar, piecewise-C 1 continuous system crosses a switch- 
ing manifold in such a way that its eigenvalues "jump" across the imaginary 
axis. This theorem requires several nondegeneracy conditions; in particular, 
the criticality of the bifurcation is determined by the linearization of the 
critical equilibrium. This is to be contrasted with the classical smooth theo- 
rem where the distinction between supercritical and subcritical bifurcations 
depends upon cubic terms. Our result extends the result given in |FPT97 
that applies to piecewise-linear systems. 

A related theorem was proved in [ZKB06J for the case of a piecewise-C 2 
system (not necessarily continuous) for the case that the equilibrium is as- 
sumed to be fixed on the switching manifold, but its left and right eigenvalues 
change with a parameter. 

2 Discontinuous Bifurcations 

Piecewise-smooth, continuous odes may contain bifurcations that do not ex- 
ist in smooth systems. For instance, if an equilibrium crosses a switching 
manifold as a system parameter is continuously varied, we expect the eigen- 
values to change discontinuously because of the discontinuity in the Jacobian. 
In this case the equilibrium may disappear or its stability may change, this 
is a discontinuous or border- collision bifurcation |Lei06j . 



2 



Consider a PWSC system that depends upon a parameter /i and suppose 
that z*(p) is an equilibrium that lies on a switching manifold at /x = 0. 
Furthermore, suppose that the switching manifold is co dimension-one and is 
smooth at the point z*(0). Then, without loss of generality, we can choose 
coordinates so that z*(0) = and the unit vector ei is the normal vector 
to the switching manifold. To reflect these choices, let z = (x,y), with x 
representing the first component, and y the remaining n — 1 components. 

As with smooth systems, knowledge of local behavior is gained by com- 
puting the Jacobian, A(z) = DF(z), of the equilibrium near the bifur- 
cation. Though the Jacobian is not defined at z*(0) = (0,0), the limits 
Al = lim x ^o_ A(x, 0) and Ar = lim x ^ 0+ A(x, 0) do exist since F is piecewise 
smooth. In this case, continuity of the system implies that all of the columns 
of Al and An must be equal except for the first. The system is PWSC if the 
first columns do indeed differ. 

Feigin gave one classification of border-collision bifurcations according to 
the spectra of Al,r [D BGIV0"2] . In particular, let sl,r be the number of nega- 
tive, real eigenvalues of Al,r, respectively. It is not hard to see that if sl + sr 
is even the equilibrium generically persists and crosses the switching mani- 
fold, while if sl + is odd the equilibrium does not cross the boundary but 
is destroyed in an analogue of the saddle-node bifurcation. Similar results are 
valid for maps where a simple criterion for a discontinuous period-doubling 
bifurcation can also be given [DBFHH99J. 

Even when the stability of an equilibrium does not change upon collision 
with a switching manifold, its basin of attraction can undergo a dramatic 
change. In the piecewise linear case, orbits of a nominally stable equilib- 
rium can be unbounded; this has been called a "dangerous border-collision 
bifurcation" [HAJN04j . 

This paper is concerned with the planar case when Ai and Ar correspond 
to focus-type equilibria of opposing stability. The discontinuous bifurcation 
arising in this situation displays similar properties to an Andronov-Hopf bi- 
furcation in a smooth system. 

As an example, consider the piecewise-linear continuous system: 

x = —x — \x\ +y , 

(1) 

y = — 3x + y — n . 

For a given value of /x, the system has a unique equilibrium, namely a stable 
focus at (— /i, — 2/x) when /x < (with eigenvalues, — | ± ^-i) and an 

unstable focus at (— 0) when /x > (with eigenvalues, \ ± ^p-i). The 
2/-axis is a switching manifold and the equilibrium crosses this manifold at 
the origin when /x = and changes stability. Fig. [T] shows phase portraits of 
([[]) for negative and positive values of /x. A stable periodic orbit exists for 
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all positive values of /i and grows in size with /i. This situation is analogous 
to a supercritical Andronov-Hopf bifurcation in a smooth system. 




3 Normal Form 

More generally, consider a two-dimensional, piecewise-C fc , continuous system 
of ordinary differential equations in z = (x,y) € M 2 , with k > 1. Assume 
that there is an equilibrium, z*(n), that crosses a switching manifold 
rameter fi is varied. If we assume this crossing occurs at a differentiable point 
on a switching manifold, and that this manifold remains locally differentiable 
under small variations of fi, we may further assume that the switching mani- 
fold coincides with, say, the ?/-axis and the crossing occurs at the origin when 
\i is, say, zero. Such a system can be written as: 

f F L (x,y;fi) , x < , , 

\ F R (x,y;v), x>0 ' ^> 

where Fl and Fr are C k vector fields with k > 1. Since the vector field 
is continuous, F L (0,y;fi) = F R (0,y; /i). We also assume that F £ (0,0;0) = 
Fr(0, 0; 0) = to ensure that z*(0) = 0. Let j e {L, i?} represent the index 
for the two pieces of (j2J), so that we can rewrite the system as z = Fj(z; fi) 
where j — L when x < and j = R when a; > 0. 

Expansion of each part of (T2]) about the origin yields 

x = p(n)fi + aj(fi)x + b(fi)y + o(x, y) , 
y = ^(A*)A* + Cj(/i)a; + d(ji)y + o(x, y) , 

where p, g, a^, b, Cj, and <i are real- valued C k ~ l functions of \i, for j e {L, i?}, 
and the correction terms go to zero faster than the first power of x or y. 
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The factor /x in each of the constant terms makes explicit the fact that the 
constant terms are zero when \x = 0. The Jacobian matrices of the left and 
right systems in ([3]) at the origin are 



Aj(ji)=DF j {p& l i) = 

We suppose A L (0) and A R (0) have eigenvalues 

A L ± iu L , and - X R ± iu R , (5) 

respectively, where \ L , \ r ,ul,u r > 0. Notice this assumption implies that 
Cj(0) and 6(0) ^ for otherwise Aj(0) would be triangular and have real 
eigenvalues. Moreover note that detA,(0) = A| + cj| > 0. Finally, since 
continuity implies that 6(/x) is the same for both matrices, the direction of 
rotation for each subsystem must be the same. 

The equilibrium of ([3]) that lies at the origin when /x = will non- 
tangentially cross the switching manifold y = if the parameter 

r=p(0)d(0)-g(0)6(0)^0. (6) 

In this case, by rescaling the parameter /x i— > T/x, we can effectively set 



6(/x) 
d{ji) 



(4) 



r = i . 



Two other simple transformations can be done to simplify the coefficients 
of Q. First, since 6(0) ^ 0, we may use the transformation y i— > 6(0)?/, to 
set 6(0) = 1. This implies that when /x is small, the rotation direction is 
clockwise. 

Second, we may perform a /x-dependent shift of y to set p(/x) = for 
small /xQ To see this, note that the implicit function theorem implies that 
there is a neighborhood of the origin in (y, /x) such that there is a unique 
C k function y(fi) that satisfies x(0,xy(/x); /x) = with y(0) = 0. This follows 
because i(0,y;/x) is x(0,0;0) = and D y x{0, 0; 0) = 6(0) = 1. Thus 
the transformation y ^ y — xy(/x), effectively sets j9(/x) = for all sufficiently 
small /x. It is easy to see that the remaining coefficients q, aj, 6, Cj, and d are 
still C k ~ x functions of /x (for small /x) and are unchanged at /x = 0. 

The system ()3]) now becomes 



x 

y 







/x + A^/x) 



X 

y 



+ o(x,y) 



(7) 



1 This simplifies the proof of Th.[T]in Sj4j in particular allowing us to center the Poincare 
map (fn?]) , at the origin, instead of at some /^-dependent point on the y-axis. 
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for small enough /i. Here we have 6(0) = 1 and since p(0) = we have 
T = — g(0) = 1 so that q(0) = — 1. Note that the eigenvalues of Al(0) and 
A R (0) are still given by §5§. 

Let x* L (fi) be the x-component of the equilibrium when it exists in the 
left half-plane, and x* R (fi) in the right half-plane. It is easy to see that 

4M = -T2^-2 + °(m) ie{L,i?} (8) 

Thus when /i is small and positive [negative] the equilibrium is located in the 
left [right] half-plane and is a repelling [attracting] focus. 



4 Andronov-Hopf-like Bifurcations 

We will show that, as with an Andronov-Hopf bifurcation in a smooth system, 
a periodic orbit of ([7]) is created at // = and grows in amplitude as /i is 
either increased or decreased, given a single nondegeneracy condition. Recall 
that the nondegeneracy condition for the smooth case is a cubic coefficient in 
the normal form |Kuz04l Wig03 . The sign of this coefficient also governs the 



stability of the limit cycle, i.e., the criticality of the Hopf bifurcation. For the 
discontinuous analogue that we are studying, the nondegeneracy condition 
depends only upon the eigenvalues of the linearized system. 

Theorem 1. Suppose that the vector field (Ej) is continuous and piecewise 
C k , k > 1, in (x,y,/j,), and has an equilibrium that transversely crosses a 
one- dimensional switching manifold when \i = at a point z* where the 
manifold is C k . Suppose further that as fi — > + the eigenvalues of the 
equilibrium approach \i ± iu^ and as fi — > 0_ they approach —Xr ± ioor, 
where \l, Xr, ujl, ojr > 0. Let 

A = — ~ — , (9) 

Ur 

denote the criticality parameter. 

Then if A < there exists an e > such that for all < /x < e there 
is an attracting periodic orbit whose radius is O(fi) away from z* , and for 
—e < fj, < there are no periodic orbits near z* . 

If, on the other hand, A > 0, there exists an e > such that for all 
—e<fi<0 there is a repelling periodic orbit whose radius is O(fi) away 
from z* , and for all < fi < e there are no periodic orbits near z* . 

In order to prove Th. [TJ we will use the transformed system (JTj) so that 
the switching manifold becomes (locally) the y-axis. The periodic orbit will 
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be obtained as a fixed point of a Poincare map, P, of the positive y-axis to 
itself. This map is obtained as the composition of two maps Pl and Pr that 
follow the flow in the left and right half-planes respectively. We must prove 
that there is a neighborhood of the origin in which these two maps are well- 
defined. A fixed point of P will be shown to exist using the implicit function 
theorem. This is essentially the same approach as that used in [MM76] to 
prove the smooth Hopf bifurcation theorem and in |FPT97j to prove the same 
result for piecewise-linear systems. 

It suffices to consider only /i > 0, for when fi < the transformation 
(x,fi,t) i — ► — (x, /i, t) produces a new system displaying the same properties 
as those listed for ([7]). The signs of fi and A become reversed and the stability 
of any periodic orbits is flipped because we are reversing the direction of time. 



Proof. Step 1: Define Pi and Pr 



As previously argued, without loss of generality we may consider the C k 
system ([7]). Consider the "left" and "right" systems (taken as though they 
were separately valid for both signs of x): 



X 







V . 




. ?0) 


X 







y 







/! + A L (fi) 
H + A r (ij,) 



x 

y 

X 

y 



+ o(x,y) 
+ o(x,y) 



(10) 
(11) 



Denote the flows of these equations by ipt(x,y;fi) and y?f(x, y;^), respec- 
tively, for (x, y, n) in some neighborhood of the origin. 

For sufficiently small d~j and <5 Mj > 0, we will define the first return map 



P R :[0,S R ] x [0,5, 



by Pn(y ; fx) = y\ where y\ is the first intersection of tpf (0, y ; jj) with the 
negative y-axis or origin for t > 0, and 



Pr 



-h,0] x [0,5, 



ML J 



by Piiyi, ji) = y 2 where y 2 is the first intersection of <pf (0, yx\ fi) with the 
positive y-axis or origin for t > 0. 



Step 2: Show P R and P L are well-defined when /i = 

When fi = the origin is a hyperbolic equilibrium of both ([TO]) and (fTTT) ; 
consequently, Pr(0,0) = Pl(0,0) = 0. Since when /i = the matrices A,- 
have eigenvalues ([5]), the trajectories of the linearized systems spiral about 
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the origin, repeatedly intersecting the y-axis. The spiral behavior is clock- 
wise because in both linear systems x(0, y; 0) = y. The same can be said of 
trajectories of ffTU]) and (TTTT) sufficiently close to the origin because, by Hart- 
man's theorem, there exist neighborhoods about the origin within which ffTOl) 
and fTTTl) are C 1 conjugate to their linearizations [Har60j . Thus when fj, — 0, 
Pr and Pl are well-defined. 

Step 3: Show P R is well-defined for small fi > 

When > 0, the origin is no longer an equilibrium of (TTTj) . The needed 
properties of Pr can be obtained by considering the local behavior of the 
trajectory that passes through the origin. This can be deduced by computing 
approximations to the flow of ( fiTl) on the y-axis. Since 

x(0,y,fu) = b(/i)y + o(y) , 

6(0) = 1, and is continuous, there is a neighborhood of the origin such 
that the positive y-axis flows into the right half-plane and the negative y-axis 
flows into the left half-plane. Moreover, since 

2/(0, y; /i) = g(/i)/i + d(fi)y + o(y) 

and since g(0) = —1 and g(/i) is continuous, then for small fj, > 0, y < 
for sufficiently small y. Thus the trajectory that passes through the origin, 
does so tangent to the y-axis and without entering the right half-plane, see 
Fig. [2J Since the flow is continuous, there exists a non-empty interval [0, Sr] 
on the positive y-axis that maps into the negative y-axis. Thus the map, 
Pr, is well-defined for some S fJiR > 0. Furthermore Pr(0; /i) = (since our 
definition of Pr allows for an intersection at t = 0). Finally Pr is C k since 
(ED is C k . 

Step 4: Show P L is well-defined for small fi > 

Since 

F L (0,0;/i) = (0,g(/i)/i) (12) 

and g(0) = —1, the initial velocity of the trajectory ^(0,0; fi) is directly 
downwards. The equilibrium of (TTUl) has the x- value given by (JSj) , and there- 
fore lies in the left half-plane when fi is small and positive. 

For small /j,, the equilibrium lies close to the origin so that we may use the 
linearization of the flow at the equilibrium to approximate (pf(0, 0; fi). Recall 
that the eigenvalues of Al(0) are ± iul, where Xl,^l > 0, thus (x* L) y* L ) 
is a repelling focus at fi = 0, and by continuity remains a repelling focus 
when /i is small enough. Consequently the flow ^(0,0; fi) initially spirals 



S 




Figure 2: Sketch of the phase portrait of © when // > 



clockwise around the equilibrium solution within the left half-plane. Since 
the equilibrium is repelling, before <pf(Q, 0; /i) has completed 360° about the 
equilibrium it will intersect the 7/-axis at some point P(0; (i) = y(fi) > 0, 
see Fig. |2j Since the flow is continuous, Pl must map a non-empty interval 
[— 5l,0] to points on the positive y-axis above Thus the map, Pl, is 

well-defined for some S^ L > 0. Also Pl is C k since f lTUl) is C*. 



Step 5: Define P and compute its derivatives at (0 + ,0 + ) 

The results above imply that there are 5, 5^ > such that the Poincare map, 
P : [0, 5} x [0, <y -> M+ defined by 

P = P L oP R , (13) 

is C fc . This map is sketched in Fig. [3j 

Note that since Pr(0; /j.) = 0, its right-sided derivative with respect to ji 
as /i — ► + is £) A4 Pr(0; 0+) = 0. Moreover, since Pl(0; /i) = and since 

y(0) = and and y(fi) = 0(/i) and is positive for /i > 0, we have 

a = DM0; 0+) = lim ^ > . 

Thus the right-sided derivative of P with respect to /i is also positive 
rvP(0; 0+) = ZVMO; 0+) = a > . 
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To compute the value of D y P(0 + ; 0), consider the linear system 



z 



a 1 
b d 



z 



with eigenvalues —A rkiu. For points on the y-axis, this system has the flow 

A ' r sin(cjt) 



¥>t(0,J/) 



— (a + A) sin(wt) + lu cos(ut) 



Consequently it maps the positive y-axis onto the negative y-axis in time 
r = - so that <p r (0,y) = (0, — e~ Xr y). Moreover, this system approximates 
the flow of (jTTI) in the right half-plane near the origin. Therefore as yo — > 0+, 

PrKVo] 0) — > — e w « ?/o- By a similar argument, as y\ — > 0_, Pt(yi;0) — >■ 

— e "i Therefore as y — ► + , 



^(y ;0) 



-e w i 



-e 



2/o 



Thus 



A,P(0 + ;0) = e " 



(14) 



Step 6: Show the periodic orbit exists when A < 



The function V : [0,5] x [0,5^] — > R defined by V(y ,fi) = P(y ;fi) — y is 
C k . Zeros of V correspond to fixed points of the Poincare map, P, and thus 
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periodic orbits in the system (JTj). In order to apply the implicit function 
theorem to V at the origin we must first smoothly extend its definition to 
a neighborhood of the origin. Let V : [—5, 5} x [—5^, <5 M ] — > K be any C fc 
function with V(yo,fi) = ^(2/o>/f) whenever yo, fi > 0. 
Note that 

i) 7(0,0) = y(0,0)=P(0;0) = 0, 

ii) D y V(0, 0) = D y P(0 + ; 0) - 1 = e An - 1, 

iii) £> M F(0, 0) = L> M P(0; 0+) = a > 0. 

Therefore, if A ^ 0, then the implicit function theorem implies that there 
is a C k function y + satisfying V(y + (fi), //) = for all /i in some neighborhood 
of fi = such that ?/+(0) = and 



; £^(0,0) l-e A - 



Consequently, when A < 0, we have D^y^O) > 0. Thus for small \i > 0, 
> so that V(y + (fi),n) = V(y+(ji),fi) = implying that y + Qi) is a 
periodic point. In other words, the graph of P(y; /i) necessarily intersects 
the diagonal as sketched in Fig. [3j and when < \i < 5^, ((7j) has a periodic 
orbit that intersects the positive y-axis at y+(fi). 

The periodic orbit intersects the negative y-axis, at, say, = Pr (?/+(//); //). 

The function y-(/i) is also C k and vanishes at // = 0. Moreover, 

D„y-(0 + ) = D y P R (0 + ; 0)Z^ + (0 + ) + D M P fl (0; 0+) 



< 



The resulting bifurcation diagram is sketched in Fig. HI The radius of the 
periodic orbits, by any sensible definition, grows linearly with respect to /i, 
to first order. 



Step 7: Show the periodic orbit is attracting 

The stability of periodic orbits can be deduced by calculating the value of 
D y P(y + (n);n). When A < 0, flUD implies that < D y P(0 + ,0) < 1. Since 
this function is C k ~ l , we therefore have < D y P(y + (fj,); //) < 1, for all 
sufficiently small \x > 0. Thus the periodic orbit is attracting, as we should 
expect since the equilibrium is repelling. 
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Figure 4: A sketch of a bifurcation diagram when A < 0. and y£(/x) correspond to 

the y-component of the equilibrium. y+(fJ.) and y-(p,) correspond to the values of the two 
intersections of the periodic orbit with the y-axis. Solid [dashed] curves denote attracting 
[repelling] solutions. 

Step 8: Show there are no periodic orbits for A > 

If A > 0, (IT51) implies that y' + (0) < 0, so that the locus of zeros, y+(fj), fails 
to enter the positive quadrant (yo, fi > 0) near (0,0). Since, by the implicit 
function theorem, y+(fi) is the unique solution that emerges from the origin, 
it follows that V(yo,n) ^ 0, for all sufficiently small yo, fi > 0. Hence in this 
case, there are no periodic orbits. □ 

Theorem [T] implies that the criticality of the discontinuous Hopf bifur- 
cation depends on sgn(A) in ([9]). It is not difficult to understand why this 
makes sense geometrically. If (JTj) has a Hopf cycle for small fi, it must en- 
circle the equilibrium and spend time in both the left and right half-planes. 
Suppose the equilibrium is repelling and lies in the left half-plane. Then, 
within the left half-plane, the Hopf cycle completes more than 180° around 
the repelling equilibrium of ( jTUl) . However, within the right half-plane it 
completes less than 180° around the inactive, attracting equilibrium of the 
right vector field (TTTT) . In order that the orbit be periodic it must, in some 
sense, spiral outward exactly as much as it spirals inwards. Consequently, 
the attracting nature of the attracting equilibrium must be stronger than the 
repelling nature of the repelling equilibrium. Since the time to move 180° 
is n/uj, this requirement is equivalent to X r /ujr > \ l /ujl, hence A < 0, in 
agreement with the theorem. 

In smooth systems, Hopf cycles grow as ellipses and the period of an ar- 
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Figure 5: A trajectory approaching a Hopf cycle for |T|) when /.i > 0. The upper 
equilibrium (blue) is the unstable focus of the left flow and the lower equilibrium (green) 
is the inactive stable focus of the right flow. The angles that the Hopf cycle subtends 
within each half-plane about each spiral center are labeled 0l and (3r. 



bitrarily small Hopf cycle is easy to calculate |Kuz04[ |Wig03| . In order to 
calculate the period for our situation, we must first understand the shape 
of the Hopf cycle. As the Hopf cycle shrinks to a point, it becomes better 
approximated by the periodic orbit in the corresponding piecewise-linear sys- 
tem. A periodic orbit in a piecewise-linear system (like (pQ)), shrinks to zero 
in a self-similar manner due to an inherent scaling symmetry. Thus an arbi- 
trarily small Hopf cycle generated by a discontinuous Hopf bifurcation takes 
the shape of the periodic orbit in the corresponding piecewise-linear system, 
hence the orbit consists of two spiral segments, see Fig. [51 Again, consider 
the case A < where the equilibrium is repelling and lies in the left half- 
plane. The Hopf cycle completes, say, (3l > n radians around the repelling 
equilibrium and 13 r < it radians around the inactive, attracting equilibrium, 
see Fig. [51 As \i — > 0+ the ratio Pl/Pr limits on a finite value strictly greater 
than one. However, even the knowledge of these angles does not give the 
period of the cycle since the time taken along each spiral segment is given 
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by two different quantities, 7l/Vl and "/r/ujr, respectively. Here we observe 
that < 7# < 7T < 7l < 2tt, but in general, (3l ^ 7l and (3r ^ jr. The 
period of the Hopf cycle is 

— (— + — + 0(fi) . (16) 

As a nonlinear example, consider the piecewise continuous system: 

x = —x — \x\+y, 

2 ( 17 ) 

y = -3x + y- n + 3y , 

which is identical to ([1]) except for the addition of a single nonlinear term. 
By theorem [TJ (I17p will display the same local dynamical behavior about the 
origin for small // as (pQ) . Fig. shows a bifurcation diagram of (ITTj) . Two 

0.2 



0.1 








-0.02 -0.01 0.01 



Figure 6: A bifurcation diagram for (|17|) . The dotted curve is the unstable focus, and the 
dashed curve a saddle. The solid curve, existing for /j, < 0, denotes a stable equilibrium. 
The solid curves for fj, > denote the two locations at which the stable periodic orbit 
intersects the y-axis. 

equilibria are born in a saddle-node bifurcation at \i — — ~ —0.02083, and 
exist for all larger values of \i. The stable node becomes a stable focus at 
[i = 2 — |v3 ~ —0.02073 and finally an unstable focus at jj, = when its 
eigenvalues jump across the imaginary axis from — | +iu>L to ~ + iujR, where 
Wi = |v3 and cj^ = Iv^ll- As predicted by Th. [TJ a stable periodic orbit 
is created at fi = exists for small fi > 0. Moreover near the origin, the 
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bifurcation diagram looks the same as that shown in Fig. HI Phase portraits 
of ( IT7I) for fi = ±0.01 are shown in Fig. El 

The periodic orbit is destroyed in a collision with the saddle equilibrium 
in a homoclinic bifurcation at fi « 0.01127. The variation of the period with 




0.004 0.008 0.012 

Figure 7: Period of the Hopf cycle of (fT7j) that is created at the origin when /i = 0. 

respect to \i is shown in Fig. [71 Note that the period at ji = is different 
from the nominal value, T(0) ~ 5.522, that it would have if 7l = 7^ = 7r in 

(USD. 




Figure 8: Phase portraits of (TTT]) when (a) [i = -0.01 and (b) [i = 0.01. 
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Co dimension- two bifurcations can arise if the genericity conditions A / 
and r / are not satisfied. As A crosses zero the bifurcation changes from 
supercritical to subcritical, and at T = the equilibrium does not cross the 
switching manifold transversely. In both cases, higher order terms are needed 
for a local analysis. For the special case that the equilibrium is fixed on the 
switching manifold, q(fj,) = 0, a periodic orbit is created if a parameter 
causes A to cross zero |ZKB06j . A further condition required for Th. [1] is 
that the equilibrium intersects the switching manifold at a point where it is 
smooth. The scenario where this is not true could perhaps be understand by 
transforming the system so that the switching manifold lies on the positive 
halves of the x and y axes; a special case where the equilibrium remains 
at a corner point was treated in [ZT05J. For this case, the parameter A is 
replaced by a sum of the ratios $l(\j)/5s(\j) of the eigenvalues multiplied by 
the opening angle of each sector. 

5 Conclusions and Open Problems 

We have shown there exist Andronov-Hopf-like bifurcations in piecewise- 
smooth continuous systems. The three major differences between this bifur- 
cation and the classical Hopf bifurcation are: 

1. An arbitrarily small Hopf cycle consists of two spiral segments as op- 
posed to being elliptical. 

2. The amplitude of the Hopf cycle grows linearly with respect to the 
system parameter instead of as the square root of the parameter value. 

3. The criticality of the bifurcation (i.e., the stability of the Hopf cycle) is 
determined by linear terms — the parameter A of (Q — instead of cubic 
terms. 

It would be nice to extend our results to the case of an n- dimensional 
system with a smooth codimension-one switching manifold such that the 
spectra of the matrices Al and Ar differ by one pair of eigenvalues. The 
difficulty here is devising a version of the center manifold reduction that is 
used in the proof of the smooth Hopf-bifurcation theorem [MM761 IKuz04] . 
The higher dimensional case is complicated by the fact that an equilibrium 
on a switching manifold can be unstable even when both Jacobian matrices 
have all of their eigenvalues in the left-half-plane; this was demonstrated for 
a three dimensional, piecewise linear example by [CFPT06] . 
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